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ABSTRACT 

We construct and evolve three one-parameter families and one two-parameter family of steady-state 
models of stellar disks embedded in live dark matter (DM) halos, in order study the dynamical and 
secular phases of bar evolution. These models are tested against those published in the literature 
in order to extend them and to include the gaseous component in the follow up paper. Specifically, 
we are interested in the angular momentum, J, redistribution in the disk-halo system during these 
two evolutionary phases without distinguishing between the resonant and non-resonant effects. We 
confirm the previous results and quantify for the first time the dual role that the DM halos play in 
the bar evolution: more centrally concentrated halos dilute the dynamical processes of the initial bar 
growth, such as the spontaneous bar instability and the vertical buckling instability, and slow down 
the J transfer, while facilitating it in the secular phase. The rate of J transfer in the disk and the 
halo is followed up in order to identify sites and times of peak activity in J emission and absorption. 
Within the corotation radius, Rcr, the disk J remains nearly constant in time, as long as Rcr stays 
within the disk — a sign that the lost angular momentum to the outer disk and the halo is being 
compensated by an influx of fresh J due to the outward motion of i?cr- We demonstrate that this is 
feasible as long as the bar slowdown dominates the loss of J inside Ra. Next, we find that in some 
models the bar pattern speed stalls for prolonged time periods, i.e., the bar exhibits a constant rate 
of tumbling when R^ is located outside the disk. This phenomenon appears concurrent with the near 
absence of J transfer between the disk and the halo, and is associated with the halo emitting J at 
the corotation resonance and absorbing it at the inner Lindblad resonance. Furthermore, we confirm 
that stellar bars generally display the corotation-to-bar size ratios in the range of 1 — 1.4, but only 
between the times of the first buckling and i?cr leaving the disk. Hence, the corotation-to-disk size 
ratio emerges as an important dynamic discriminator between various stages of barred disk evolution. 
Finally, we analyze a number of correlations between the basic parameters of a barred disk and a halo, 
some already reported in the literature and some new. 

Subject headings: galaxies: evolution - galaxies: halos - galaxies: kinematics and dynamics - galaxies: 
spiral — galaxies: structure - stellar dynamics — cosmology: dark matter 



1. INTRODUCTION 

Bars are expected to have a profound effect on disk 
galaxy evolution because they constitute a major depar- 
ture from an axial symmetry and hence facilitate angular 
momentum and mass redistribution. Despite prolonged 
and focused investigation of bar formation and evolution 
many theoretical and observational issues remain unre- 
solved. Among these are the origin of bars in the universe 
— spontaneous or tidally induced, evolution of bar frac- 
tions with redshifts, and theoretical predictions of bar 
pattern speeds. The intricacies of the bar growth or de- 
cay with time are obscure — how do bars capture orbits? 
Do all bars extend to their corotation radii? Lastly, vari- 
ous aspects of nested bar systems are only now beginning 
to be analyzed. 

Because the underlying dynamics of bars is strongly 
nonlinear, their numerical modeling has spearheaded the 
efforts to understand their orbital structure, ability to 
channel the angular momentum to the outer disk and 
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dark matter (DM) halo, efficiency in triggering the ra- 
dial gas inflows, and more. While we have gained some 
insight to all of these processes, much remains to be done 
and one expects additional effects to surface. 

In this work we aim at deeper understanding of an- 
gular momentum redistribution in the disk-halo system 
facilitated by the stellar bars in the presence of the gas 
component. Here, we analyze a set of equilibrium colli- 
sionless models which differ by a single parameter from 
each other. In the associated paper, we add the gas com- 
ponent to the system without affecting the mass distribu- 
tion there (Villa- Vargas et al.. Paper II, in preparation). 
Our coUisionless modeling is in a way complementary to 
that of Athanassoula (2003). We reproduce some of the 
correlations between the bar, disk and halo properties 
discussed in the above work and quantify additional cor- 
relations found. All these will be compared to models 
with gas. 

Importance of a galactic spheroidal component in bar 
evolution has been gradually understood, from its sup- 
posedly stabilizing effect on the disk (Ostriker & Peebles 
1973; Efstathiou et al. 1982), to serving as disk angu- 
lar momentum sink (e.g., Sellwood 1980; Weinberg 1985; 
Debattista & Sellwood 2000; Athanassoula 2002), and 
appears to be full of controversies (e.g., Christodoulou, 
Shlosman & Tohline 1995a,b; Athanassoula 2008). Study 
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of the dominant role of resonance interactions between 
the bar and surrounding orbits was pioneered by Lynden- 
Bell & Kalnajs (1972) and applied to bar and halo orbits 
by Tremaine & Weinberg (1984), Athanassoula (2002), 
Martinez- Valpuesta, Shlosman & Heller (2006), Cevcrino 
& Klypin (2007), Weinberg & Katz (2007a,b) and Du- 
binski, Berentzen & Shlosman (2009). Toomre (1981) 
argued that the bar growth can be damped by the in- 
troduction of an inner Lindblad resonance (ILR), by 
cutting off the swing amplification mechanism for the 
TO = 2 modes. Hotter stellar disks are expected to be 
more stable as well (e.g., Athanassoula & Sellwood 1986). 
Counter-intuitively, in the long run, a disk embedded in 
a more centrally concentrated halo develops a stronger 
and longer bar (Athanassoula & Misiriotis 2002). 

A comprehensive analysis by Athanassoula (2003) has 
shown that the efficiency of resonances in angular mo- 
mentum transfer depends on the DM and stellar disper- 
sion velocities and the DM densities in the vicinity of 
the major resonances. As the resonances sweep across 
the phase-space, about 20% — 30% of DM particles are 
locked in the lower order resonances at any given time 
(Dubinski et al. 2009). Additional aspects of the bar 
evolution are known to be influenced by the spheroidal 
components, such as the vertical buckling instability in 
the bar (Berentzen et al. 2007). 

In general, J transfer in the disk-halo systems can 
depend on three parameters, namely, the particle pop- 
ulation near low level resonances and the velocity dis- 
persions in the disk and halo (e.g., Athanassoula 2003). 
In principle, this means dependence on the halo central 
mass concentration and on radial velocity dispersions in 
the disk and the halo. We, therefore, focus on how pa- 
rameters of the DM halo distribution, such as its mass, 
concentration, extent, and rms velocities affect the disk- 
halo evolution and the angular momentum transfer medi- 
ated by the stellar bar. Specifically, we are interested in 
how the presence of gas in the disk influence the angular 
momentum redistribution in the system. For the pur- 
pose of clarity, we separate our discussion of coUisionless 
models (this Paper I) from models with gas (Paper II). 
CoUisionless models are used in order to verify the basic 
details of bar evolution in order to serve as benchmarks 
for models with gas. They also are used to calibrate our 
results to those of Athanassoula (2003). 

Our paper is structured as follows: §2 describes the 
numerics, initial conditions and model parameters. §3 
provides results on basic evolution of stellar bars in our 
sets of models and analyzes the redistribution of angular 
momentum in these systems, while §4 is focused on corre- 
lations between various parameters in the bar-disk-halo 
models. 
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2.1. Initial conditions 



The initial conditions were created following the pre- 
scriptions and density profiles from Hernquist (1990). 
The mass volume density distribution in the disk is given 
in cylindrical coordinates by 



Pd{R,z) 



exp{~R/h) sech^ 



20 



(1) 



where is the disk mass, /i is a radial scale length and 
zq is a vertical scaleheight. 
The density of the spherical halo is given by 
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where Af^ is the mass of the halo, is a Gaussian cutoff 
radius and 7 is the core radius, a is the normalization 
constant defined by 



{l-^/^gexp(g2)[l_erf(g)]}-l 



(3) 



with q = 7/rc. 

The particle velocities, dispersion velocities and asym- 
metric drift corrections have been calculated using mo- 
ments of the coUisionless Boltzmann equation. Since 
models thus constructed are not in exact virial equilib- 
rium, the halo component was relaxed for t '-^ 40 in the 
frozen disk potential. 

2.2. Model parameters 

Our goal is to investigate the effect of the DM mass 
distribution on the evolution of stellar bars in a live disk- 
halo system. We choose a specific model, defined here as 
the standard model (SD), where the DM density profiles 
are specified by three parameters, Mh, 7 and r^, corre- 
sponding to the DM halo mass, DM core radius and the 
DM Gaussian cutoff radius. Based on the SD model, 
three sequences of models have been created, each se- 
quence resulting from varying one of the parameters only. 
In the fourth sequence, we simultaneously modify two pa- 
rameters at the time, 7 and rc, in attempt to target the 
outer halo only, while trying to keep the inner halo un- 
changed. In the SD model, the halo-to-disk mass ratio 
in the inner R — 0.6 is kept to unity. Table [T] lists the 
values of the SD parameters. 



TABLE 1 

Parameters of the standard model 



2. NUMERICS AND MODELING 

We use the TV-body part of the FTM-4.4 hybrid code 
(e.g., Heller & Shlosman 1994; Heller, Shlosman & 
Athanassoula 2007b) to evolve the stellar disks and DM 
halos. The gravitational forces are calculated using the 
FalcON routine (Dehnen 2002) which scales as 0{N). 
The units of mass and distance are taken as lO^^M© 
and 10 kpc respectively. This makes the unit of time 
as 4.7 X 10^ yr when G = 1, and the velocity unit 
208 km s^^. The gravitational softening is egrav = 0.016 
for stars and DM particles. The models consist of a stel- 



HALO 
Parameter Value 


DISK 
Parameter Value 


A^DM 




Af. 


2 X lO'"^ 




3.15 




0.63 


n 


8.55 


Rt 


1.71 


7 


0.1425 


h 


0.285 




2.85 


20 


0.057 






Q 


1.5 



Note. — Q is the Toomre parameter fixed at i? = 2.4h, where 
h is the thickness of the disk; rt and Rt are numerical truncation 
radii in the halo and the disk. All values are given in dimensionless 
units, §2. 
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The initial mass volume density profile in the disk is 
kept unchanged in all models. The spatial distribution of 
the stellar particles is thus identical in all realizations of 
the initial conditions. On the other hand, the velocities 
of the stellar particles have been adjusted in each models 
to provide for the changing rotational support against 
the combined gravitational potentials of the varying disk- 
halo systems. The four model sequences are as following: 

Mh sequence: the mass of the halo, Mh, has been re- 
duced to 70% and 40% (models M70 and M40 re- 
spectively) of the SD model. 

7 sequence: the DM core radius 7 was increased by a 
factor of 2 and 4 with respect to the SD model 
(models C30 and C57 respectively). 

rc sequence: the DM Gaussian cuttoff radius Vc was de- 
creased to 77% and 56% of that in the SD model 
(models T22 and T16 respectively). 

-Mh — Tc sequence: in these hybrid models, MTl and 
MT2, both Tc and Afh have been varied, to keep 
the inner r < 0.1 DM density profile close to that 
of the SD model. 



TABLE 2 

Halo parameters for the model sequences 



Model 7 rc 

"SD 3.15 0.1425 2.85 

M70 2.20 0.1425 2.85 
M40 1.26 0.1425 2.85 

C30 3.15 0.30 2.85 
C57 3.15 0.57 2.85 

T22 3.15 0.1425 2.20 
T16 3.15 0.1425 1.60 

MTl 2.40 0.1425 2.20 
MT2 1.71 0.1425 1.60 
Note. — All values are given in dimensionless units, §2 

The DM halo parameters in each model are listed in 
Table[21 The number of particles and the numerical trun- 
cation radii are the same as in the SD model. The halo 
mass Mh is conserved in sequences 7 and rc. The DM 
density profiles of all models after relaxation in the frozen 
disk potential are shown in Figure [H The resulting DM 
density profiles exhibit monotonic decrease along the Afh 
sequence, show a progressively larger core in the 7 se- 
quence, move the outer halo mass inward to the inner 
halo, across r ^ 1.4, in the sequence, and increase 
the outer halo mass while leaving the core unchanged in 
the hybrid models sequence. The initial circular veloc- 
ity curves, showing the contributions from the disk, the 
halo and the total, are displayed in Figure [2] Models 
M40, C30 and C57 host maximal disks which dominate 
the potential of the inner part. Models SD, M70, MTl 
and MT2 have equal contributions of disk and halo at 
t = 0, while T22 and T16 are halo dominated. With 
the development of the stellar bar, the disk typically be- 
comes even more dominating, although the inner halo is 




-1 
R 



Fig. 1. — The DM halo density profiles at t = after relaxation in 
the frozen disk potential. Each panel shows the model sequences. 
From top to bottom: Afh, 7, f'c and the hybrid Afh — rc sequences. 



also dragged in an adiabatic contraction (e.g., Dubinski 
et al. 2009). 

3. RESULTS 

The model sequences described in Table 2 differ only 
by the parameters of DM distribution as given by Eq. 2. 
These include the total mass of the halo, its core size and 
the Gaussian truncation radii. As the bar properties and 
angular momentum evolution are both heavily depen- 
dent on the mass distribution in the system, the choice 
of the free parameters allows us to fine-tune the changes. 
Because the DM appears to be on the receptive side of 
the angular momentum transfer, we target the halo as a 
whole, and its inner and outer parts separately. Discus- 
sion and comparison with published models is made in 
§4. 

3.1. Bar strength and pattern speed evolution 

To gauge the strength of the bar we use the Fourier am- 
plitude A2 of the TO = 2 mode normalized by the m = 
mode. It is obtained by various methods — here we show 
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Fig. 2. — Circular velocities of all models at t = after halo 
relaxation in the frozen disk potential. Each panel displays the 
stellar disk rotation curve (red, long dashed), halo (blue, dashed) 
and the total curve (black, solid). 



the results of integration over all the disk, and integra- 
tion over restricted cylindrical volumes of a particular 
interest. We define A2 in the integration limit over the 
fixed radial range R ~ 0.1 — 1.8. Second, we define 
when integrate over R = 0.1 — range, where Rh is the 
bar size defined in §3.2. Intuitively, the second definition 
reflects the bar properties more fully because is not 
diluted by the disk properties. Both A2 and A2b show 
some similarities, e.g., the same peaks, raises and drops, 
and differences, e.g., the relative strength of the peaks 
(Fig-ISl first and second columns). 

Most models show common stages in the evolution of 
A2 and A2b. (1) an initial exponential steep rise, (2) a 
subsequent peak or plateau followed by a sudden drop, 
ensued by (3) a more gradual and sustained rise, ending 
by (4) a saturation of the bar strength. The duration of 
each stage varies from model to model, but stages "1" 
and "2" are much shorter than stage "3". We refer to 
stages "1" and "2" as the dynamical evolution of the 
bar, and stages "3" and "4" as the secular evolution of 
the bar. Models M40 and T16 are the exceptions. The 
former model has the least massive halo and exhibits a 
prolonged secular growth of the bar after the first peak. 
The bar amplitude does not reach saturation even in the 
Hubble time. The latter model has the most massive halo 
core and the bar strength never shows a clear first peak 
seen in all other models. However a careful analysis of A2 
and especially A2h reveals that the exponential growth 
is terminated at t ~ 120 and the subsequent evolution 
lucks the vertical buckling instability (see §3.4), but oth- 



erwise follows the path of other models. We take a closer 
look at these trends in §4. We note that our decision to 
create one-parameter sequences naturally leads to a wide 
spectrum of bar properties and can include some of the 
extreme behavior. 

Some clear trends of the bar strength behavior can 
be observed along each of the sequences defined in §2.2. 
First, in all models, except the two exclusions mentioned 
above, A2 and A2h have saturated by t ^ 230. Because 
A2b is tailored for the bar, the saturation of this param- 
eter is more obvious (more about this in §3.3). Second, 
making the inner halo less massive (i.e., less centrally 
concentrated), either by reducing the total halo mass or 
by increasing the size of the central core (sequences Mh 
and 7, respectively), results in a shorter rise time of the 
bar instability and hence brings up the bar earlier. Third, 
increase in the mass of the inner halo and decrease in 
the outer one, has the most dramatic effect on the bar 
strength, substantially increasing the timescale of the bar 
instability. This hints at DM mass concentration being 
important rather than the total halo mass. Varying the 
mass of the outer halo alone has a much smaller effect 
on the bar. Finally, in all one-parameter sequences, the 
first A2 or A2h peak forms a progressively more extended 
plateau for less concentrated halos, respectively. Again 
T16 is an exclusion. 

Furthermore, A2 and A2h show high frequency vari- 
ability in some models (e.g., Figure [J]), being strongest 
in models SD, C30 and T22. This variability can also be 
seen in i?b and Q^, evolution, the bar length and its pat- 
tern speed. This variability is limited to the time of the 
bar secular growth and dies out with saturation of the 
bar strength, A2 and A2h- These oscillations coincide in 
time with the presence of four short arms located close 
to the bar end, and two spirals in the outer disk beyond 
the bar radius. Both sets of spirals have pattern speeds 
larger than that of the bar. The presence of more than 
one strong pattern of waves or modes in the disk can 
be the result of a non-linear mode couphng (e.g.. Tagger 
et al. 1987; Sygnet et al. 1988) that gives rise to beat 
waves of a frequency Wbcat = wi ± 0^2 and the azimuthal 
wave number mbeat = "nii ± m2, where mi and m2 are 
the (azimuthal) wave numbers of the interacting waves 
and LOi and L02 are their frequencies, to — mil, where fi's 
are their pattern speeds. This type of coupling has been 
observed and quantified in numerical simulations (e.g., 
Masset & Tagger 1997; Martinez- Valpuesta 2006). De- 
battista & Sellwood (2000) has also observed the high 
frequency variability, but no quantitative analysis of the 
mode coupling was attempted. 

We test whether the bar and the outer two spirals, 
both TO = 2 modes, couple to give rise to the observed 
beat mode TObeat — 4 arms. We verify that the frequency 
fj^bcat of the TOboat = 4 mode is actually the sum of the 
frequencies of the interacting modes. Table [3] shows the 
frequencies t^b = 2r2b, Ws ~ 2r2s and Wboat obtained 
with the mode coupling analysis, with an error of ±0.02. 
Here fls is the pattern speed of the outer to = 2 spirals. 
The measured frequency of the beat mode (column 5) 
matches the frequency predicted by the non-linear cou- 
pling (column 6). We can explain the oscillations in A2 
as resulting from superposition of the major axis of the 
bar with the outer spirals, which are the second strongest 
mode (after the bar itself) present in the disk. The fre- 
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Fig. 3. — Time evolution of normalized bar amplitudes, A2 (left column) and (middle column), and the bar pattern speed ^2^ (right 
column) in all models. The model sequences are indicated in the upper left corners. All the data has been smoothed with a high frequency 
Fourier filter. 




Fig. 4. — Close up of the A2 evolution in the SD model, showing 
the high frequency oscillations during i ~ 70 — 160. No smoothing 
has been applied to the data. 



quency /a2 of the A2 oscillations must thus be equal 
to the frequency at which the spirals align with the bar 
/sb = (ij-'s — Wb)/27r. As can be seen from columns 7 and 
8, these two frequencies coincide within the error limits. 

Figure [3] shows the evolution of fib in our models, 
which is well defined by the time A2 ~ 0.05. In all 
models, the initial value of Sib is very similar, with the 
scatter of ±0.1. Although the details of the evolution of 



TABLE 3 
Mode coupling in the disk 



Model 


Time 








+ L^B 


/a2 


/sb 


SD 


106 


1.80 


4.86 


6.66 


6.66 


0.49 


0.49 


C30 


113 


1.32 


4.24 


5.58 


5.56 


0.45 


0.46 


T22 


149 


1.80 


5.72 


7.58 


7.52 


0.61 


0.62 



Note. — (1) the model; (2) the time at which the measurements 
were taken; (3) oj^,; (4) lJb', (5) ttibeat! (6) the expected frequency of 
the beat mode ui^, + lUb', (7) measured frequency of A2 oscillations 
(8) frequency of alignment of bar and spirals. 



r^b differ from model to model, most importantly, often 
shows an (anti)-correlation with A2 and , as empha- 
sized already by Athanassoula (2003). This is not always 
observed in cosmological simulations of disk evolution 
(Heller, Shlosman & Athanassosula 2007a; Romano-Diaz 
et al. 2008). We find, however, that the anti-correlation 
between the bar strength and its pattern speed are lim- 
ited to times prior to saturation of A2 and A2h- In §3.3 
we show that this corresponds to times when the bar 
corotation radius lies within the disk. In three models, 
M70, T22 and MTl, the A2-ilb correlation is maintained 
even after the corotation leaves the disk — in these mod- 
els the total J of the disk is conserved during this time 
period of flat A2 and fib- 

A close examination of Sib evolution in Fig. [3] reveals 
prolonged periods of Slb^ const. These are found either 
in the dynamical phase, prior to the first buckling but 
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when the bar is already sufficiently strong (SD, M70, 
C57) , or during the late stage of secular evolution, espe- 
cially pronounced in M70, T22 and MTl, and to a lesser 
degree in C57. These latter cases are of most interest to 
us. A similar behavior was discussed before by Valen- 
zuela & Klypin (2003), who related it to an abnormally 
low dynamical friction of the bar against the background. 
We return to this issue in §3.6. 

3.2. Evolution of the bar length 

The semi- major axis of the bar, i?b, is defined here 
as the radius where the bar equatorial ellipticity drops 
by 15% off its peak. Martinez- Valpuesta et al. (2006) 
has tested this method by contrasting it with the most 
reliable way of determining the bar size using the last 
stable orbit supporting the bar (Martinez- Valpuesta et 
al. 2006). It seems robust when applied after the first 
maximum of the bar strength. The ellipticity of the bar 
at different radii is obtained by fitting ellipses to the 
isodensity contours in the face-on disk. Figure [5] displays 
the evolution of i?b for all the models. 

Similarly to the evolution of A2 and A2h, the bar length 
has initial period of a fast growth that reaches maximum, 
sometimes followed by a drop and subsequent sustained 
but slower rise, in most models. Some models show a 
degree of saturation, toward the end. The first peak in 
i?b always coincides in time with that of A2, except in the 
model T16 which lacks it. The peak in i?b, however, is 
not as pronounced as that in A2 and there are additional 
caveats. As in A2, the initial growth varies with a model 
in such a way that, in models with the less centrally 
concentrated halos, the time span of the initial growth is 
shorter. Furthermore, in some cases, additional drop in 
the bar length can be seen after some time of the secular 
growth, in the time range of t ~ 150 — 260. This drop 
appears to be related to the formation and disappearance 
of the ansae at each tip of the bar (Martinez- Valpuesta 
et al. 2006). The ansae eventually 'detach' from the 
bar and becomes misaligned with its major axis. The 
transition to this misalignment results in a fast decrease 
of the bar length, taking place at t ~ 180 for model SD, 
- 190 for M70, 22 for C30, 190 for C57, 205 for T22, and 
165 for MTl and MT2. 

In general, evolution in the bar length does not neces- 
sarily go in tandem with the changes in A2 and is much 
less monotonic. For example, the drop in i?b around 
t ~ 180 in the SD model has no counterpart in the A2 or 
A2h evolution. Of course, a stronger A2 does not mean 
a longer bar, as seen in the last At ^ 80 of the evolution 
of the sequence Afh models. 

Overall, the A/h and Tc sequences show a substantial 
dispersion in the bar sizes at the end of the simula- 
tions, while two other sequences end up with very similar 
bars. In all sequences, the bars differ substantially during 
the intermediate times, especially after the first buckling 
time. 

3.3. Evolution of the bar corotation radius 

The bar corotation radius i?cr has been computed us- 
ing linear approximation (the second column in Fig. [5]). 
It grows with time as a consequence of the bar slowdown. 
i?cr keeps growing throughout the simulation, the excep- 
tions are models M70, T22 and MTl, in which a plateau 
is reached. We use this characteristic radius in order 



to construct and follow the evolution of two important 
ratios, namely, Rcr/Rh and Ra/Rd- 

The position of i?cr has a profound effect on the angu- 
lar momentum transfer from the bar region to the outer 
disk and the DM halo (Athanassoula 2002; Martinez- 
Valpuesta et al. 2006). The DM particles at the coro- 
tation resonance are responsible for much of the angular 
momentum absorption by the halo. Hence, the avail- 
ability of halo or outer disk resonant particles is of the 
prime importance to the evolution of the bar. Specifi- 
cally, it is important to know whether i?cr lies within the 
stellar disk in our models at all times. We define Rd as 
the disk radius which encloses 98% of the disk mass, and 
the evolution of the ratio Rcr/Rd is shown in Figure [5l 
At the time of the bar formation, Ra /Rd lies between 
0.5 and 0.8 in all models. This ratio raises with time 
enough to move the corotation out of the disk with the 
exception of models M40 and T16, where bars grow at 
the slowest pace during secular evolution. The growth of 
the bar seems to be sensitive to the moment at which i?cr 
reaches the disk edge — even if the growth continues, it 
proceeds at a slower pace. This change happens with a 
slight but measurable delay of At ^ 20. Moreover, as 
noted before, the detachment of the ansae from the bar 
happens just before Rcr/Rd ^ 1, which manifests itself 
by a sharp drop in the bar length. 

We pay a special attention to the bar evolution when 
Rci/Rd ^ 1- Fig. [S] displays an interesting behavior of 
this ratio in the above regime — it flattens with time 
before resuming its growth. Analyzing this behavior, we 
find that as the i?cr moves outside the disk, some stellar 
particles become trapped at i?cr- m = 2 spiral arms 
appear connecting the outer edge of the disk to these 
trapped particles. This is characteristic of all models 
where the above conditions exist. At t ^ 200 — 250, 
depending on the model, the spirals dissolve, and their 
particles form an amorphous cloud outside the disk. We 
return to this issue in §3.6. 

The bar must lie of course inside its i?cr, because orbits 
outside i?cr do not support the bar. The morphology of 
the offset (gas) shocks in the numerical bars has been 
argued to constrain the ratio Rcr/Rh, restricting it to 
Rci/Rh = 1-2 ± 0.2 (Athanassoula 1992; see also early 
work by van Albada & Sanders 1982). Limiting our dis- 
cussion of this ratio only to times when Rcr/Rd ^ 1, we 
observe it dropping to the range ^ 1.0—1.4 after the first 
buckling, and staying generally below ^1.8 until the i?cr 
moves out of the disk. Model T16 (no buckling!) does 
not follow this trend and rather shows a monotonic decay 
to ~ 1.9 at the end of the run. On the other hand, M40 
— the only other model with Rcr always in the disk, dis- 
plays a remarkably constant Rcr/Rh = 1-2 ± 0.2 almost 
over the entire run (Figure [5]). Naturally, after i?cr moves 
out of the disk, this ratio is bumped to above 2. In Pa- 
per II, this ratio will be tested against the shape of the 
offset shocks in the gas. 

3.4. Vertical buckling in the bar 

The A^-body bars undergo phases of vertical buck- 
ling instability where transient asymmetries develop in 
the form of vertical oscillations of the equatorial plane 
and asymmetric thickening of the disk (e.g.. Combes & 
Sanders 1981; Combes et al. 1990; Pfenniger & Friedh 
1991; Raha et al. 1991), and results in the formation 




of the so-called peanut /boxy-shaped bulges (Kormendy 
& Kennicutt 2004 and refs. therein). These phenom- 
ena can be recurrent and affect different parts of the bar 
(Martinez- Valpuesta et al. 2006; Berentzen et al. 2007). 
We measure the buckling asymmetry with the Fourier 
coefficient Ai^^ of the m — 1 component in the rz-plane, 
again normalized by the m=0 mode. In all our mod- 
els, except TI6, we detect at least a single phase of the 
vertical buckling. The buckling always happens simulta- 
neously with the first A2 drop (e.g., Martinez- Valpuesta 
& Shlosman 2004), which we confirm here. In all the 
model sequences, models with less centrally concentrated 
halo exhibit earlier and stronger buckling. This happens 
because A2 reaches its first maximum earlier in these 
models. We do not find any clear correlations between 
the strength and time of the second buckling and the 
initial properties of the models. In all cases of repeated 
bucklings, i?cr lies within the disk. 

We note that T16 does not exhibit the drop in A2 
or A2h, but does show a break in their slopes around 
t ^ 120. Moreover, while all our models develop 
peanut/boxy bulges abruptly after these drops, T16 de- 
velops this bulge gradually over the secular stage of its 
evolution. This model appears as a nice example of a 
secular buildup of peanut/boxy bulges as a result of the 
diffusive action of the vertical inner Lindblad resonance 
(Friedh & Pfcnniger 1990; Martinez- Valpuesta & Shlos- 
man 2004). Berentzen et al. (2007) have shown that 
mass concentrations resulting from gas accretion to the 
center damp the vertical buckling and alter the bulge 
shape, making it more elliptical, if the gas fraction in 
the disk is high. T16 is a purely collisionless model with 
the most centrally concentrated halo. It is possible that 
a qualitatively different behavior is the corollary of this 
concentration. Clearly, the action of the vertical reso- 
nance leads to a secular and not dynamical buildup of a 
peanut-shaped bulge in this case. 



We also note that in all models the nonlinear inner 
Lindblad resonance (ILR) first appears only after the first 
buckling, although the linear analysis claims their exis- 
tence from the start (Pfenniger & Friedli 1991; Martinez- 
Valpuesta et al. 2006). The appearance of this resonance 
is the consequence of the increase in the central mass 
concentration, by a factor of ~ 2, as a result of the bar 
buckling, which drags inward the DM as well (Dubinski 
et al. 2009). 

3.5. Evolution of bar-to-disk mass ratio 

The bar increases its mass by capturing particles from 
the disk within its i?cr- In Figure [7] we plot the fraction 
of the stellar mass contained inside the bar, AfbarA^d- 
To estimate the bar, we sum the particle masses within 
a rectangular box aligned with the bar major axis, with 
dimensions given by the major and minor axes of the bar 
and \z\ < 0.5. 

By the end of the run, the bar has captured between 
75% and 85% of the disk mass (the exception is model 
T16 whose bar growth has been severely delayed by the 
halo mass concentration). The initial stage of bar for- 
mation is associated with an intense capture of disk par- 
ticles. Eventually this rate declines, and even though 
the bar mass keeps growing it does so at a much more 
modest pace. The growth of Afbar/Afd clearly correlates 
with the evolution of the bar size in Fig. [5] (left frame) — 
an increase / decrease in one of them corresponds with an 
increase/decrease in the other. Typically, the growth of 
the bar size and of Mbar/-Md saturates when the Rcr/Rd 
is driven above unity. However, surprisingly, there are 
exclusions of this behavior, e.g., C57. 

We do not find correlation between the bar mass and 
A2. It may happen that the bar keeps capturing mass 
while A2 is saturated, as in model SD at i ^ 210 and later 
on. In summary, the ratio Mx^^T-jM^ saturates around 0.8 
in in all models, except T16. 
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Fig. 6. — Evolution of the bar vertical asymmetry given by 
the Fourier coefficient of the m = 1 mode in the r^-plane 

corotating with the bar major axis. No filtering was applied to this 
data. 



3.6. Evolution of the angular momentum 

The backbone of the bar is composed from material 
moving along eccentric orbits of relatively low angular 
momentum, compared to the circular orbits of the same 
energy. Moreover growing bars consume more circular 
orbits, i.e., those with larger J/E ratio, redistributing 
their angular momentum to the outer disk and halo. The 
capacity of different regions of the disk and halo to emit 
and absorb angular momentum is thus of a prime im- 
portance to the bar evolution, and largely determine its 
ability to grow by trapping additional disk orbits. Rather 
than estimate the effect of the resonances on this process 
(Athanassoula 2002; Martinez- Valpuesta et al. 2006), we 
focus on the total J transfer within the disk and to the 
halo, i.e., accounting for the resonant and non-resonant 
J redistribution. In Figure [5] we plot the evolution of 
the z component of the total angular momentum of the 
disk, Jd, as well as in the cylindrical volumes centered on 
the symmetry axis of the disk and separated by i?cr, i-c, 
R <i?cr ('^d,in) and i?cr £ R (>/d.out), and similarly for 
the halo (Fig. [5]). Since the total angular momentum in 
each model is conserved, the evolution of Jh in the halo 
is an inverted mirror image of that in the disk. 

The initial angular momentum, Jd(i = 0), of the disk 



Fig. 7. — Evolution of bar-to-disk mass ratio M\^^j- / M^. 



varies from model to model due to the differences in the 
halo mass distribution. It changes very little before the 
vertical buckling in the bar, while Jd.in and Jd,out an- 
ticorrelate. This means that the halo gains a negligible 
amount of J, which flows nearly exclusively to the outer 
disk during this phase (see also Fig. [5]). After the buck- 
ling, Jd and Jd,out drop monotonically till the end of the 
run. This is not true about Jd,in which stays approxi- 
mately constant until the time when i?cr leaves the disk 
(Fig. [5J see also Martinez- Valpuesta et al. 2006). How 
this is possible? 

While it is not surprising that the loss of J within the 
corotation leads to an outward motion of i?cr, the appar- 
ent 'conspiracy' in Jd,in~ const, over such an extended 
time period is somewhat puzzling. It hints at some reg- 
ulation mechanism which fine tunes the loss of angular 
momentum to the outer disk and to the halo by gravita- 
tional torques with its influx across the corotation. We 
attempt to answer this question, using obvious simplifi- 
cations: the explanation of this phenomenon lies in the 
intricacies of Jd,in balance and the resulting behavior of 
rib- The angular momentum within i?cr is 

Jd,i„ = 27r/ dRR^j:{R,t)vt, (4) 
Jo 

where vt = oiVc is the average tangential velocity of par- 
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Fig. 9. — Evolution of the z components of the halo angular 
momenta, Jjj (solid, black) and Jjj ^ut (dashed, green), enclosed 
by the cylindrical volumes R < Ra and R > Ra, and z = ±2. 



tides inside Ra, which also defines the coefBcient a. We 
require Jd.in^ const, between the times of buckling^ and 
when Rc-r/Rd ^ 1, and take the time derivative of Jd.in 
using the Leibniz formula, 



1 dJds, 
2tt dt 



R^-[j:{R,t)vt]dR 



+v,Rl^{R„,t) 



dRc, 
dt 



= 0, 



(5) 



where E is the disk surface density which can be obtained 
from eq. [T]by integrating over z. We assume Vc ~ const. 
This is justified because contribution of J within the ve- 
locity turnover radius is small and at larger radii, Vc is 
dominated by the isothermal halo and is nearly indepen- 
dent of radius (Fig. [2]). The first term in eq. [5] describes 
the emission of J by the inner disk, within an instanta- 
neous Rci, due to the torque T imposed by the DM and 
the outer disk. For brevity, we assume that the stellar 
bar extends to i?cr and this torque acts upon it. In other 
words, all the material within i?cr resides in the bar. The 
second term in eq. [5] represents the influx of Jd,in due to 
the advance of Rcr and the resulting addition of a new 
mass within this radius on nearly circular orbits. Next, 
we relate the motion of i?cr to the slowdown of the bar, 
dRcr/dt = Rcr = —Rci-flh/^b- The torque on the bar 
gives the rate of change of its angular momentum (and, 
therefore, change in the angular momentum of the inner 

^ In fact, we can impose this condition already from t = 0, 
as variation of is relatively small during the buckling (e.g.. 

Fig. [9)1 
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disk): 



J ni 



(6) 



where lb is the moment of inertia of the bar and Jcirc is 
the angular momentum of the internal circulation within 
the bar. (The second term does not appear in eq. 10 of 

Athanassoula (2003).) Hence, T — tl^Ih + ^^b^b + >^circ- 
Assuming that the dominant loss of the angular momen- 
tum by the bar is due to the slowdown of its tumbling 
f2b (not a trivial assumption and definitely not a general 
one!), we can re- write eq. [5l replacing its first term by 
the action of T, as 



"b 



RlAh) = 0, (7) 



where we have taken Mb ~ i?^j.S]cr- If the second term 
in eq. [7] loosely represents the moment of inertia of the 
bar, the net change in Jd,in is indeed negligible, explain- 
ing its near constancy in Fig. [8] (mid column) up to the 
time when i?cr leaves the disk. In summary, we show 
that the constancy of the angular momentum within _Rcr 
can be indeed explained if a number of straightforward 
assumptions is made to allow for an analytical estimates. 
These assumptions are as following: (1) the bar extends 
to near i?cr during this time and the disk mass within 
this radius lies mainly in the bar; (2) the rate of change 
of the angular momentum in the bar due to the change in 
its moment of inertia and internal circulation is smaller 
than that resulting from the bar slowdown; (3) the circu- 
lar velocity in the disk around the region and beyond 
is independent of radius; and (4) the average tangential 
velocity of the material in the bar is a fixed fraction of 
the disk circular velocity. While this behavior shows up 
in all our models, as long as _Rcr lies within the disk, a 
more general set of mass distributions should be tested 
before making a final conclusion. 

While the halo absorbs all J emitted by the disk, it 
is instructive to look into which parts of the halo are 
especially active in this process. Figure O displays the 
evolution of the absorbed Jh^n and Jh.out- We observe 
that the inner halo becomes more receptive to the process 
with time, except in M40 and T16. This is clearly related 
to the outward motion of i?cr- We note also, that the rate 
of J absorption by the inner halo increases substantially 
after Ra /Rd ^ 1, and the halo serves as the only sink of 
the angular momentum. Asymptotically, one can divide 
the models in two classes: when the rate of J absorption 
by the inner halo goes to zero or even becomes negative 
(e.g., M70, T22, MTl), or when it declines only shghtly 
(e.g., SD, M40, C30). 

An alternative way to detail the exchange of the an- 
gular momentum in the disk-halo system is to divide 
the disk and halo into a number of concentric cylindrical 
shells. We have constructed a 2-D map of the angular 
momentum in each shell as a function of radius and time. 
It is more revealing, however, to plot the time derivative 
of Jd and Jh at each radius. This analysis was performed 
separately for the halo and the disk. In Figures fTOl to fT4l 
we display color-coded diagrams of the rate of change of 
angular momentum Jd,R = (dJd/dt)^ with radius and 
time for each run, and the same for the halo. Red color 
corresponds to emission and blue to absorption of the 
angular momentum. All the diagrams corresponding to 
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Fig. 10. — The rate of angular momentum flow J as a function of 
radius and time for the SD model (see also the text). The top row 
corresponds to the halo and the bottom row to the disk. Red/blue 
colors represent the absorption/emission of J using a linear scale 
in color. The main resonances are indicated by the solid lines: the 
outer/inner Lindblad resonances (OLR/ILR), the outer/inner ul- 
traharmonic resonances (OUHR/IUHR) and the corotation (CR). 
The cylindrical shells have AR = 0.1 and z = ±2. Time smearing 
is At = 3. 



disks have been calibrated identically. DM halos have 
also been calibrated uniformly among themselves. 
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Fig. 11. — J as a function of radius and time for models M70 
(left column) and M40 (right column). The top row corresponds 
to the hal o a nd the bottom row to the disk. Color code and lines 
as in Fig. [lO] 

The most striking features in Figures [10] to [14] is the 
continuity of color, i.e., emission and absorption bands, 
both in the disk and the halo, and their outward shift. 
The bands appear nonuniform in time and space. The 
halos are almost exclusively absorbers of the angular 
momentum, while the disks show regions of absorption 
and emission. From the halo evolution in the SD model 
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Fig. 13. — J as a function of radius and time for models T22 
(left column) and T16 (rigt column). The top row corresponds to 
the ha lo a nd the bottom row to the disk. Color code and lines as 
in Fig. [TO] 



(Fig. [TT|) . we note that Jh is absorbed preferentially at 
three distinct bands which shift outwards with time. 
Other and weaker absorbing regions are observed be- 
tween them and at larger radii. The absorption/emission 
bands in the disk are alternating, unlike in the halo. 

To relate the emission/absorption bands to dynamical 
characteristics of the disk-halo system, we specify the po- 
sitions of linear resonances on the color map. These are 
used for guidance purpose only, as the positions of linear 
resonances maybe misleading, especially for the strong 
bars, and especially for the inner Lindblad resonance(s). 

We find a close correspondence between the behav- 
ior of Jd(R) bands and the location of the linear reso- 
nances (superposed solid lines), namely, the corotation. 



inner/outer Lindblad resonances (ILR/OLR), and the in- 
ner/outer ultra- harmonic resonances (UHR). While it is 
clear that the redistribution of J in the disk and the halo 
is a strongly nonlinear process, some general conclusions 
can be made. 

First, as expected, the emission and absorption be- 
tween the inner disk and the inner halo components 
strongly correlate, with a certain interference of the outer 
disk fFigs. [TTI - fH)) . Here we refer to inner/outer with re- 
spect to the Rci-. The peak emission of J by the disk 
happens at the band associated with the ILR. A weaker 
emission band lies at the inner UHR (see also Fig. 10 
of Martinez- Valpuesta et al. 2006). The strongest ab- 
sorption in the disk lies at the corotation and extends 
to the outer UHR. Weaker absorption is associated with 
the OLR in the disk. The halo absorption bands are cen- 
tered on the ILR, inner UHR/corotation and the outer 
UHR/OLR. However, it is nonnegligible in other regions 
as well. We note that producing horizontal slices in the 
halo reveals that the lower halos, i.e., Az — ±0.3, absorb 
mostly at the ILR while the upper halo absorb mostly 
at Rcr- This is a refinement to previous results which 
shown that the main absorption is at the CR resonance 
(Athanassoula 2002; Martinez- Valpuesta et al. 2006). 

Second, we observe temporal correlations between the 
emission and absorption peaks in the disk and the halo. 
Over the length of the run, both the disk and the halo 
experience from two to three 'times of activity,' except 
T160. The first such activity time is clearly associated 
with the first peak in A2 and the subsequent buckling. 
The other activity time lies toward the saturation of A2 , 
late in the evolution of the bar. Some models, like SD, 
show an additional activity in between. This seems to be 
associated with i?cr crossing i?d — the halo absorption 
peaks at the corotation during this time. The resonances 
move outwards faster during intense J absorption. The 
outward migration of the resonances is of course related 
to the bar slowdown, resulting in the resonances sweeping 
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across tlie_ disk and halo material. 

Third, Jd(R) and Jii(R) in Figs. fTTUTH appear to corre- 
late with Jd,in, ^d,out, Jhdn and Jh.out in Figs. [5] and [H It 
also correlates with A2 and especially with A2h- Lastly, 
the maxima of disk emission correspond to the maxima 
in the halo absorption. 

Next, we return to the atypical behavior of some stel- 
lar bars in their secular phase, which exhibit f2b~ const., 
discussed in §3.1. In all the models where we observe this 
phenomenon, the ratio Rcr/ Rd > 1 at the time of detec- 
tion. This behavior, as expected, is associated with the 
rate of / transfer from the disk to the halo: in all these 
models, J —^ 0, both in the halo absorption or disk emis- 
sion during fib ~ const. (Figs. [TT] to [H]). The outward 
motion of the disk resonances stalls as well (Figs. [5] and 
[TT] — [T4|) . Moreover, there is an indication that this is 
concurrent with the reversal of J flow: while the disk 
still shows a weak emission at its ILR, the halo shows 
a weak emission at i?cr and a strong absorption at the 
ILR. The latter remains the most active resonance in the 
halo. The emission of J by the halo at _Rcr is rather un- 
usual and has not been observed before in any models in 
the literature, to the best of our knowledge. 

Lastly, the disk inspection reveals the presence of 
two trailing spiral arms extending to i?cn and when 
Rcr/Rd > 1, connecting the outer edge with _Rcr (see 
§3.3). Some stellar particles appear to be trapped at i?cr 
in its outward motion up to t ~ 250, depending on the 
model. These particles appear to concentrate at specific 
azimuthal locations of the CR "circle" with respect to 
the bar (in the bar frame of reference). 

The particles in the spirals are not trapped. These 
spirals are prominent in Figs. [TTl [T4l and are responsible 
for emission at i?cr and absorption of J in the lUHR-CR 
band, as i?cr moves outside the disk. While the num- 
ber of these particles is not large, their specific angular 
momenta is the largest in the disk. 

4. DISCUSSION: BAR-DISK-HALO CORRELATIONS 

We have constructed and evolved three one-parameter 
and one two-parameter sequences of DM halos hosting a 
standard asymmetric stellar disk. All models developed 
bars and we have followed their evolution over a Hubble 
time. The properties of the bars, such as their strength, 
pattern speed, size, corotation radius, ratios of corotation 
to bar lengths, and more, have been related to those 
of the disks and host halos. Based on our results, we 
now test a number of correlations found by Athanassoula 
(2003). We also obtain additional ones — all of these 
will be used in order to compare with the gas models 
(Paper II). 

For this purpose, we deal with three entities: the in- 
ner, bar- forming disk (defined within i?cr), the disk as a 
whole, and the DM halo. We adopt the halo-to-disk mass 
ratio, Mh/Md within the disk velocity turnover radius, 
-Rturn = 2.2/1 (e.g., Sackett 1997) at i = 0, as an inde- 
pendent variable, where h is the disk radial scalelength 
defined in section \2A\ In the following, we also use the 
radial and vertical dispersion velocities in the disk and 
radial dispersions in the halo, and test how the efficiency 
of J redistribution within the system depends on these 
parameters. 

To quantify the role of the halo in the bar instability, 
we calculated the exponential timescale of A2 growth. 
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Fig. 15. — Average rate of angular momentum loss, (| J|), by the 
inner stellar disk during its dynamical (upper) and secular (lower) 
phases of the bar evolution as a function of the M^^/M^ ratio. The 
rate (J) has been calculated within -Rturn (upper) and within i?d 
(lower). The mass ratio has been calculated within -Rtum • 

The 

averaging was performed over the time period from t = to the 
end of an exponential growth of the stellar bars (upper) and from 
the time of the first minimum of A2, i.e., the of the buckling, to 
its saturation, t = 213 (about 10 Gyr) (lower). The green circles 
represent the two hybrid models. The attempted linear fit in the 
lower frame does not include the lower right corner point (T16). 



^ban against M\JMd within i?turrn which represents the 
halo concentration, normalized by the disk mass within 
the same radius. A linear correlation was obtained, 
namely, that an increase in the DM mass concentration 
leads to a slower bar growth before the first buckling. 
However, the bar evolution in the secular phase cannot 
be described by an exponential growth. Looking for a 
more universal measure of the bar growth at all times, we 
decided in favor of angular momentum change in the in- 
ner and full disks. Figure [15] displays the average rate of 
angular momentum change, (| J|), in the dynamical and 
secular phases of the bar evolution, within -Rturn and i?d 
respectively. The former choice of a fixed radius (i?turn) 
results from total Jd being nearly unchanged before the 
buckling — J is exchanged predominantly between the 
inner and outer disks across i?cr at this stage. We have 
refrained from using i?cr as an inner/outer disk separa- 
tor here because its outward motion brings in fresh, high 
J/E material which 'contaminates' the J transfer due to 
the resonant and non-resonant interactions between the 
inner and outer disks. The latter choice of i?d in the 
secular phase results from Jd,in being nearly constant — 
J is mostly exchanged between the disk and the halo. 
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Fig. 16. — A2 at t ~ 213, corresponding to 10 Gyr, as a function 
of the angular momentum change in the whole disk over At = 213. 
The green circles represent the two hybrid models. 



Figure [T51 shows that clear correlations of (|J|) with 
Mh/Md persist over the Hubble time — more concen- 
trated halos provide an increasingly efficient damping of 
the bar instability by reducing the rate of J transfer to 
the outer disk (upper frame). They do facilitate such a 
transfer in the secular phase (lower frame). Significantly, 
these correlations extend down at least to Mh/Md ~ 0.5 
— the disk-dominated models. We note that T16 ex- 
hibits a qualitatively and quantitatively different behav- 
ior and was not included in the attempted linear fit to 
the secular evolution (lower frame). We have discussed 
this model in §3.4. 

The explanation for such a dual role of a DM halo in 
the bar evolution lies in determining what serves as a sink 
of Jd.in from the inner disk. As shown in §3.6, the outer 
disk beyond i?cr absorbs nearly all of J during the bar in- 
stability. A more concentrated halo during this phase will 
dilute the disk gravity and, therefore, will act against the 
resonant orbit coupling (i.e., resonant torques) between 
the inner and outer disks. On the other hand, during the 
secular phase of the bar instability, after the first buck- 
ling, the halo serves as the sink of J from the disk. In this 
case, if increase of the DM mass density leads to a con- 
current increase in the DM phase space density near the 
resonances, it facilitates the bar-halo resonance coupling, 
unless counterbalanced by a 'hotter' halo. The (J) rates 
appear similar during the dynamical and secular phases. 
But the duration of each phase can differ considerably. 
Therefore, one expects that the amount of J acquired by 
the halo will be much more substantial during the secu- 
lar phase compared to the angular momentum lost by the 
inner disk (within a fixed radius!) during the dynamical 
phase. 

It is known that the DM halo affects the bar evolution. 
In particular, it was shown that the DM halo concentra- 
tion anti-correlates with the bar growth time. Athanas- 
soula (2003) has also demonstrated the (anti)-correlation 
between fib and A2. Our Figure [T5l quantifies the dual 
role played by the DM halo during the bar evolution and 
ties this explicitly to the rate of the angular momentum 
transfer within the disk-halo dynamical system. 
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Fig. 17. — Saturated A2 sd, t = 213 (~ 10 Gyr), as a function 
of the halo-to-disk mass ratio within Rturn at t = 0. Two models 
do not reach the saturation and are shown as lower limits: M40 
(left) and T16 (right). The green circles represent the two hybrid 
models. 
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Fig. 18. — Dependence of Ai^z amplitude on the ratio of vertical- 
to-radial velocity dispersions, (ctz d/^'r.d)^- T16 is excluded as it 
does not buckle. The green circles represent the two hybrid models. 
These velocity dispersions have been measured in a ring defined by 
0.1 < R< 0.2 and \z\ < 0.1. 



Stronger bars are expected to be more efficient in redis- 
tribution of the angular momentum in the disk-halo sys- 
tems. In our collisionless models, bars strengthen mono- 
tonically only prior to the first peak and after the first 
buckling. Despite the overall non-monotonic behavior of 
the bar, its final strength correlates with the amount of 
AJd lost by the disk and acquired by the halo (Fig. [TB)) . 
in agreement with Athanassoula (2003). 

For the same reason, we expect and obtain correlation 
between A2 at t ^ 213 (-- 10 Gyr), and Mh/Md ratio 
within i?turn, i-e., the halo mass concentration (Fig. [T7)). 
thus confirming Athanassoula (2003) claim. By this 
time, nearly all modeled bar strengths reach saturation 
(e.g.. Fig. [3]), except in two models, M40 and T16, shown 
as lower limits. As shown in §3.3, A2 and saturate 
shortly after Ra/Rd ^ 1- These models represent two 
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Fig. 19. — Secular evolution of the bars: dependence of (A2), 
the average growth rate of the bar, between tbuckl a-nd t = 213 
(10 Gyr), vs cr^^d/'^c, the ratio of the radial dispersion in the disk 
to the circular velocity at i?turn. The latter ratio is calculated at 
*buckl- T16 is omitted as it does not buckle. The green circles 
represent the two hybrid models. 
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Fig. 20. — Dependence of Ai^z peak on the central halo density 
at t = 0. The green circles represent the two hybrid models. 



extreme trends in the halo mass concentration. T16 has 
the most concentrated halo, does not buckle and its secu- 
lar evolution differs from other bars. M40 represents the 
least concentrated halo with a fractionally most massive 
disk (Figs. [T] and [5]). It has the shortest exponential 
growth for the bar, ibar ^ 10- Therefore, it is puzzling 
why its bar strength does not saturate. While M40 ex- 
hibits buckling and subsequent drop in A2 and A2h, its 
secular growth is the slowest of all models and compara- 
ble to that of T16. Its bar is the longest of all models as 
the ansae do not disappear. Furthermore, its i?cr lies al- 
ways within the disk, while i?cr/^b~ 1-2 and remains flat 
for a Hubble time. All these properties appear unique in 
our models. 

While the behavior of T16 is partially understandable, 
that of M40 requires more explanation. It is difficult 
to explain within the context of the halo mass concen- 
tration only, and seems to require additional parameter. 



We invoke the disk radial dispersion velocity, o-R.d, which 
in a way also governs the bar vertical buckling instabil- 
ity through its ratio to the vertical velocity dispersion, 
Cz.d- We test this idea by first plotting the amplitude 
of buckling, Ai^z, vs this ratio calculated just before the 
first buckling (Fig. [T8|). Indeed, M40 has the smallest 
(cz.d/cR^d)^ among our models, and, therefore, is ex- 
pected and, in fact, displays the largest buckling ampli- 
tude Ai^z (Fig. H]). A side effect of this buckling is the 
overall heating of the disk. Consequently, the disk will 
be less susceptible to the bar growth in the secular phase. 
Fig.IIHshows the growthrate of the bar strength averaged 
over the secular phase, {A2), i.e., between the buckhng 
time, tbucki, and the saturation time, t 213 ('^ 10 Gyr), 
vs CTR^d/wc measured at tbucki (i-e., at the minimum of 
A2). A clear correlation exists between these parameters 
— 'hotter' disk impairs the bar growth in the secular 
phase of its evolution. Note that M40 has the hottest 
disk in this stage, which would explain its anomalously 
slow bar growth. The same correlation against CTR^d looks 
weaker and shows much larger dispersion of individual 
models. 

As mentioned before, in each separate sequence, the 
Ai z peak is stronger for models progressively less dom- 
inated by the halo. It is interesting that an anticorrela- 
tion between the initial central density of the halo and 
the strength of the first Ai^^ peak holds for our mod- 
els altogether, regardless of which halo parameter is var- 
ied (Figure [201) • Moreover, the anticorrelation exists be- 
tween the strength of the first Ai^^ peak and the ratio 
of the vertical-to-radial dispersion velocities, shown in 
Figure [m 

To summarize, we confirm the previous results and 
quantify for the first time the dual role that the DM 
halos play in stellar bar evolution (Fig. 15): more cen- 
trally concentrated halos slowdown dynamical processes 
in the disk, such as spontaneous bar instability and ver- 
tical buckling instability, as well as angular momentum 
redistribution in the system. They reverse this trend and 
facilitate the angular momentum transfer during the sec- 
ular stage of bar evolution, following the buckling. We 
follow the angular transfer in the disk-halo systems by 
varying one basic parameter at the time in order to iden- 
tify the sites and times of intense angular momentum 
flows. While we confirm the earlier works which have 
identified the ILR and corotation resonances as being 
primarily responsible for J emission by the disk and ab- 
sorption by the halo, we also find few caveats. 

First, the total angular momentum in the disk is largely 
conserved before the first buckling, thus J flows across 
the corotation to the outer disk, although the amount 
is relatively low compared to the subsequent exchange. 
Second, during the secular stage of the bar evolution, 
the angular momentum within the corotation resonance 
is largely conserved: the loss of J by this region due 
to the gravitational torques on the bar (resonant and 
non-resonant) is compensated by the influx of new ma- 
terial, rich in J, across the corotation, due to its out- 
ward motion. We also elaborate under what conditions 
this 'conspiracy' law operates. In view of this, the ra- 
tio of Rcr/ Rd emerges as important dynamic discrimina- 
tor between various paths in barred disk evolution — its 
value can be determined if bar pattern speeds are known. 
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Third, wc find that in some models the bar pattern speed 
stalls becoming nearly constant for prolonged time peri- 
ods. All stellar bars which show this behavior have their 
corotation lying outside the disk. The disk-halo angu- 
lar momentum exchange nearly vanishes and the halos 
display a weak emission of J at the corotation and ab- 
sorption at the ILR. While otherwise the bar pattern 
speeds are generally strongly decaying over the secular 
timescale, the addition of the gas component can modify 
this trend and soften the bar braking. We postpone our 
conclusions on this issue to Paper II. We also find that 
Rci/Rh ratio stays within 1 — 1.4 range, only occasionally 
spiking above it to ~ 1.8. This behavior is typical after 
the first buckling and as long as R„/Ra ^ 1- 

Lastly, we confirm some known correlations (Figs. 16 
and 17) between the basic parameters of the disk-halo 
system, e.g., between the final bar strength and MijAI^ 
ratio. Model M40 with the least concentrated halo does 
not follow this trend and we test the possibility that an 
additional parameter plays the role in this, i.e., the dis- 
persion velocities ratio crz^d/fR.d- Because M40 experi- 
ences the strongest biickling, the disk is heated up, which 
impairs the subsequent bar growth and can explain the 
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behavior of M40 in the previous correlation. Figs. 18, 19 
and 20 display new correlations between various param- 
eters. The bar average growthrate over its secular evolu- 
tion time decreases with increasing dispersion velocities 
in the disk, VE.,d/vc- We also show that in a closely as- 
sociated relation — the angular momentum lost by the 
disk, AJd, over the evolution correlates with the final 
bar strength. Finally, the amplitude of the first buckling 
depends on the central density in the DM halo. These 
correlations will be followed up when the gas component 
is present. We expect substantial modifications in the re- 
lations between various basic parameters in the disk-halo 
system. 

We are grateful to Ingo Berentzen for comments on 
the manuscript, and to Ron Buta for answering our ques- 
tions. We also thank Ingo Berentzen and Inma Martinez- 
Valpuesta for help with numerical issues. This research 
has been partially supported by NASA/LTSA/ATP, 
NSF and STScI grants to I.S. STScI is operated by the 
AURA, Inc., imder NASA contract NAS 5-26555. I.S. is 
grateful to the JILA Fellows for support. 



Hernquist, L. 1993, ApJ, 356, 359 

Kormendy, J. & Kennicutt, R.C. 2004, ARAA, 42, 603 
Lynden-Bell, D. & Kalnajs, A.J. 1972, MNRAS, 157, LI 
Martinez- Valpuesta, I. & Shlosman, I. 2004, ApJ, 613, L29 
Martinez- Valpuesta, I. 2006, Ph.D. Thesis, University of 
Hertsfordshire 

Martinez- Valpuesta, I., Shlosman, I. & Heller, C.H. 2006, ApJ, 637, 
214 

Masset, F. & Tagger, M. 1997, A&A, 322, 422 
Ostriker, J.P. & Peebles, P.J.E. 1973, ApJ, 186, 467 
Pfenniger, D. & Friedli, D. 1991, A&A, 252, 75 
Raha, N., Sellwood, J. A., James, R.A. & Kahn, F.D. 1991, Nature, 
352, 411 

Romano-Diaz, E., Shlosman, I., Heller, C.H. & Hoffman, Y. 2008, 

ApJ, 687, L13 
Sackett, P.D. 1997, ApJ, 483, 103 
Sellwood, J.A. 1980, A&A, 89, 296 

Toomre A. 1981, The Structure and Evolution of Normal Galaxies, 

CUP: Cambridge, Fall S. M. & Lynden-Bell D. (eds.). Ill 
Tremaine, S. & Weinberg, M.D. 1984, MNRAS, 209, 729 
Valenzuela, O. & Klypin, A. 2003, MNRAS, 345, 406 
van Albada, T.S. & Sanders, R.H. 1982, MNRAS, 201, 303 
Weinberg, M.D. 1985, MNRAS, 213, 451 
Weinberg, M.D. & Katz, 2007a, MNRAS, 375, 425 
Weinberg, M.D. & Katz, 2007b, MNRAS, 375, 460 



